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"^ , Abstract 

A combination of recent numerical and theoretical advances are applied to ana- 
lyze the scaling behaviour of the site-diluted Ising model in two dimensions, paying 
special attention to the implications for multiplicative logarithmic corrections. The 
analysis focuses primarily on the odd sector of the model (i.e., that associated with 
^ magnetic exponents), and in particular on its Lee- Yang zeros, which are determined 

to high accuracy. Scaling relations are used to connect to the even (thermal) sector, 
and a first analysis of the density of zeros yields information on the specific heat 
and its corrections. The analysis is fully supportive of the strong scaling hypothesis 
and of the scaling relations for logarithmic corrections. 



1 Introduction 

The Ising model in two dimensions with uncorrelated, quenched random-site or random- 
bond disorder is a paradigm for the study of the statistical mechanics of disordered sys- 
tems. With moderate dilution, the random system exhibits a phase transition different to 
that of the pure system, and the nature of this transition has been investigated for over 
three decades. 

For the disordered Ising models in two dimensions, two scenarios have arisen. The 
strong universality hypothesis maintains that the leading critical exponents remain the 
same as in the pure case and that the disorder induces multiplicative logarithmic correc- 
tions to scaling, while the weak universality hypothesis favours dilution-dependent leading 
critical exponents. While the former is now mostly favoured, especially in the random- 
bond Ising model (RBIM), the debate regarding the scaling behaviour of the specific 
heat has persisted, especially in the site-diluted Ising model (RSIM), focusing on whether 
this is characterised by an extremely weak double-logarithmic divergence or a finite cusp. 
Indeed, according to the Harris criterion pQ, the vanishing of the specific-heat critical 
exponent a in the pure Ising model, marks the borderline between the a > case, where 
disorder is relevant and where the critical exponents may change as random quenched dis- 
order is added, and the a < scenario, where this type of disorder does not alter critical 
behaviour and the critical exponents are unchanged [T]. In this borderline circumstance, 
logarithmic corrections to scaling (distinct from the logarithmic divergence in the specific 
heat in the pure Ising system) may arise. 

The issue of scaling in the RSIM is addressed here, as that is the more contentious 
version of the model. All of the thermodynamic information concerning a statistical 
mechanical system is contained in the locus and density of its partition function zeros [2] . 
In Ref . [3] , the first numerical analysis of the Lee- Yang zeros in a disordered system (the 
RSIM) was performed and the leading critical behaviour of the lowest lying Lee- Yang zeros 
was extracted through a finite-size scaling (FSS) analysis. A robust numerical technique 
to construct the density of zeros from simulational data was given in Ref. jlj. Finally, 
a self-consistent scaling theory which links the exponents characterising the logarithmic 
corrections to scaling was recently presented in Refs. [5J E|. 

Here, these numerical [3], Hj and analytical [6] advances are combined to unambigu- 
ously determine the leading scaling behaviour and potential multiplicative logarithmic 
corrections in the specific heat through the density of Lee- Yang zeros and the scaling 
relations. This and all multiplicative logarithmic correction-to-scaling exponents for FSS 
of odd functions are zero. I.e., there are no multiplicative logarithmic corrections for the 
FSS of the magnetic susceptibility, the individual Lee- Yang zeros and the density of zeros. 
This comes about through the delicate manner in which the exponents of the logarithms, 
which are nonzero in thermal scaling, balance each other out. In this way, it is established 
that the Lee- Yang zeros of disordered systems can be precisely determined numerically, 
a density-of-zeros analysis is applicable to such a system, also at the level of logarithms, 
new scaling relations for the logarithmic corrections [5], [6] in this model are confirmed 
and a negative exponent for the contentious specific heat or its multiplicative logarithmic 
correction is made unlikely. 



2 Logarithmic Corrections and Scaling Scenarios 

Recently a self-consistent scaling theory for logarithmic-correction exponents has been 
presented [HI E]. Denoting the reduced temperature by t and the reduced external field 
by h, this theory deals with the circumstances where, in the absence of field, the specific 
heat, magnetization, susceptibility and correlation length scale respectively as 

Coo(t) ~ |£Hln|£|r, (2.1) 

moo(t) ~ l^llnltll^, (2.2) 

Xoo(*) ~ \tr\\n\t\\\ (2.3) 

£«,(*) ~ |tnin|t|r. (2.4) 

At t = the magnetization scales with reduced field as 

rrioo{h) ~hi\hi.\t\f. (2.5) 

The locus of Lee- Yang zeros along the imaginary h-axis is parameterized by r and the 
Yang-Lee edge, which marks the end of their distribution, is denoted by rYx(t). In the 
symmetric phase, this scales as 

r YL (i)~|t| A |ln|t|| A , (2.6) 

while the density of these zeros along their locus (r > ryi) at criticality (t = 0) behaves 
as 0] 

9oo(r) ~r° 2_1 |lnr| a2 . (2.7) 

The leading critical exponents are related by the standard scaling relations for contin- 
uous phase transitions [7]. The theory presented in [HI E] relates the exponents of the 
logarithmic corrections in an analogous manner. 

With the strong universality hypothesis, the leading critical exponents for the dilute 
Ising models are identical to their pure counterparts: 

« = 0, 0=1, 7 = ^, 5 = 15, v = l, (2.8) 

with gap exponent 

A = /3 + 7 = y. (2.9) 

The latter equation has been verified in Ref. [3j. The exponent a^ characterizing the 
leading behaviour of the density of Lee- Yang zeros in Eq. (12.71) is [5 J 

a 2 = ^, (2.10) 

which gives a 2 = 16/15 if the strong hypothesis (J2.8P and (12.91) holds. 

Shalaev [8j and later Shankar and Ludwig [9] , and then Jug and Shalaev [10] , used field 
theory, bosonization techniques and conformal invariance to derive theoretical predictions 



for the logarithmic-correction exponents in the random-bond case with sufficiently small 
quenched dilution. These corrections, which we term the SSLJ exponents, are [5J [91 [10] 



a = 0, B=~. 7 = -, 5 = 0, = - (2.11) 

' H 16 ' ' 8 2 v ; 

with specific heat diverging as a double logarithm [H [HI [9] [11] 

Cooit) ~ |ln| In |t||| . (2.12) 

The questions addressed in the literature over the past three decades concerned (i) the 
validity of the strong scaling hypothesis and the veracity of the leading exponents (12. 8p 
in the diluted cases, (ii) the validity of the theoretical derivation of the SSLJ correction 
exponents (|2.11j) . which relies on certain assumptions regarding the nature of the dilution 
in the random-bond case (see, also, Ref. p2] for the random-site case), (iii) whether these 
exponents are dilution independent, (iv) the vanishing of a and the validity of the double 
logarithm in (I2.12p . and (v) whether or not these sets of exponents also apply to the 
random-site version of the model. 

For the two-dimensional random Ising models, an alternative scenario to (i) has per- 
sisted in the literature. This is the weak universality hypothesis and maintains that 
certain leading critical exponents change continuously as the concentration of impurity 
defects is increased [13]. In particular, the exponents a, (3, 7 and v are maintained to be 
dilution dependent while 5, rj and the ratios (3/v and j/u remain independent of the di- 
lution. Furthermore, while agreeing with the SSLJ double-logarithmic form (12.121) for the 
specific heat, Dotsenko and Dotsenko (DD) used the renormaliztion group to predict [TTJ 

Xoo(t) ~ t- 2 exp (-c(ln (- hit)) 2 ) , (2.13) 

where c is a constant related to the concentration of disorder. While the SSLJ and DD 
predictions for the susceptibility differ substantially, it is fair to say that after much 
work by various authors the strong hypothesis is now mostly favoured (especially in the 
random-bond model). 

The validity of the logarithmic-correction exponents (12.111) proved harder to establish 
quantitatively. The first direct, clear, quantitative validation of the SSLJ prediction for 
the magnetic susceptibility (that 7 = 7/8) came in Ref. [3] through series expansions. 
While the detailed scaling behaviour of (12.131) has long been ruled out, the weak versus 
strong controversy persisted. In contrast to the susceptibility, DD and SSLJ agree on 
the double logarithmic behaviour of the specific heat. However, this has been notoriously 
difficult to confirm numerically, and has been the source of much controversy. 

Simulational works, generally supportive of the vanishing of a and a and the specific 
heat diverging as a double logarithm viz. (12.121) . are found in Refs. [TS] HSl 113 Q2J and 
Refs. [T81 EH ED] for the bond-disordered and random-site Ising models, respectively (see 

also Refs. [I21[ai2DI221l231l21ll251l2Sll2Z!)- Indeed > P lots of the measured specific heat 
as a function of the double logarithm of the lattice extent are contained in Refs. [T5l [18] 
and Refs. (IBj [!§] EOj for the RBIM and the RSIM, respectively. However, in [28] it was 
claimed that such apparent double-logarithmic FSS behaviour does not necessarily imply 



Table 1: Selection of recent works supportive of the weak or strong scaling hypothesis. 





RBIM 


RSIM 


Support strong universality hypothesis 
and theoretical support for a = a = 
or numerical support for a = a = 


[11 M [221 [231 M 
|Sl El HDl E] 

dsi esi Hung 


J23I2SII2Z! 

[12] 

[THEOT 


Support for weak universality hypothesis 
and theoretical support for finite C^it) 
or numerical support for finite C^t) 


[29J 
[311 [32] 

[28, S3] 


[30] 

[33J 

[35, 36, 37, 38J 



divergence of the specific heat, and numerically based counter-claims that the specific 
heat remains finite (so that a < or a = and a < 0) in the random-bond [2S1 El] and 
random-site models [53 ESI E3 EE] also exist (see also Refs. [221 ESI ED E21 E3] ) . The 
situation is summarized in Table 1. 

The difficulties in unambiguously discriminating between the weak and strong scenar- 
ios on the basis of finite-size data were highlighted in Refs. [2HE71E9]. Indeed, in Ref. [H] 
a fit to 

Coo(t) oc \\nt\ & , (2.14) 

was attempted, and it was observed that a decreases from the value of 1 (which cor- 
responds to the pure model) as the strength of disorder is increased. This could be 
interpreted as supporting almost any reasonable value a < 1. In Ref. [37] it was pointed 
out that specific heat data in the literature, which were stated to be supportive of ( 12.121) . 
can often equally be fitted to (12.1 ft with negative a, still consistent with the Rushbrooke 
relation. 

In this paper we address the problem from a fresh perspective, namely that of partition 
function zeros (see also Ref. [3])- In Ref. [5] the scaling relation 



A = /3-7, 

for the logarithmic correction to the FSS of the Lee- Yang zeros was derived, 
diluted Ising models in two dimensions, this leads to the prediction 



A = 



15 
16 



-0.9375 



(2.15) 



For the 



(2.16) 



and verification of this value for the logarithmic correction exponent characterising the 

scaling of the Yang-Lee edge in ( 12. 6p is one of the aims of this work. 

In addition to testing the applicability and efficacy of the Lee- Yang- zero technique in 

random models, also at the level of logarithmic corrections, a central aim of this paper is 

to measure the density of zeros. Indeed, the correction exponent for the density of zeros 

is given in Ref. [5] as 

7 A + A 7 
a 2 = . (2.17) 



From the scaling relations for logarithmic corrections, in the special circumstances which 
prevail the diluted Ising models in two dimensions given in Ref. [6], S2 is related to the 
specific-heat correction exponent a appearing in (12. ip via 

A 16 * , 

a = l + 2- + a 2 = l + — A + a 2 , (2.18) 

A 15 

having used the established value (I2.9P for the leading gap exponent [3]. The elusive 
specific heat scaling exponents a, a can thus be measured from the density (12.71) together 
with (I2.10p and (12.181) and, in this way, one can distinguish between the competing a = 0, 
a = and a < or a < scenarios. 

Since it is more contentious, we address the site-diluted version. Unlike the self-dual 
random-bond version, the location of the critical temperature in the site-diluted Ising 
model is not exactly known and has to be estimated numerically. In this work, the highly 
accurate measurements for the critical temperatures reported in Ref. [19] for different 
values of the site dilution are used. 



3 Simulation of the RSIM 

The partition function for a given realization of the RSIM in a reduced magnetic field h 
is 

Z L (/3, h) = 53 exp /5 5Z e i e j a i a j + h J2 e i a i > (3- 1 ) 

Wi} V (ij) i J 

where L denotes the linear extent of the lattice and the sum over configurations {a{\ 
is taken over Ising spins Oi G {±1} and where 6j are independent quenched random 
variables which take the values unity with probability p and zero with probability 1 — p. 
For simulational purposes a regular (square) lattice with periodic boundary conditions is 
used. The percolation threshold for such a lattice in the thermodynamic limit occurs at 
p = p c = 0.592 746 . . ., so that for p < p c the lattice fragments into finite-size systems on 
which no true transition can occur [40J. Writing 

S = Y^ titjViVj , M = Y ti^i > (3-2) 

(ij) i 

and 

p L ((3;M)=J2PL(S,M)ex V (3S, (3.3) 

s 

where the spectral density pi{S, M) gives the relative weight of configurations with given 
values of S and M, the partition function in imaginary field ih is 

Z L ((3, h)=Y, Pl(/3; M) exp (ihM) = Z L (J3, 0) (cos (hM) + % sin (KM)) , (3.4) 

M 

where the expectation value has real measure. Assuming the Lee- Yang theorem holds 



Table 2: The number of samples simulated for each lattice size L at each value of the 
site-occupation probability and corresponding (3 C value (from Ref. [19]). 



L 


32 


48 


64 


96 


128 


196 


256 


p = 0.88889 p c = 0.53781 
p = 0.75 p c = 0.77125 
p = 0.66661 /9 C = 1.10 


1000 
1000 
1000 


1000 
1000 
1000 


1000 
1000 
1000 


600 

1000 

920 


600 

1000 

600 


250 
550 
270 


250 
270 
160 



[2j 13], since odd moments of the magnetization vanish for t > (/? < /3 C 
given realization of disorder are given by the values of h for which 

(cos hM) = . 



the zeros for a 



(3.5) 



Finally, for each value of L and p, these zeros are averaged over realizations of disorder 
and the resulting j th Lee- Yang zero is denoted by Tj(L). Errors associated with the zeros 
are computed as sample-to-sample fluctuations. 

We have simulated, using the Wolff single-cluster algorithm |41| . three different values 
of the dilution, namely p = 0.88889, 0.75 and 0.66661. From Ref. [19] the values of 
the critical temperatures for these three dilutions are (3 C = 0.53781(2) for p = 0.88889, 
P c = 0.77125(8) for p = 0.75, while p c = 1.10 corresponds to p — 0.66661(3). In each of 
these three cases we have run lattices of extent L = 32, 48, 64, 96, 128, 196 and 256. The 
number of samples simulated for each lattice size and each value of the site-occupation 
probability (and corresponding /3 C value) is given in Table [2j 

We have monitored the behavior of the non-local observables (such as the susceptibil- 
ity) with the Monte Carlo time by using a standard logarithmic binning of the dynamical 
data. We have checked that, in all the cases, the non-local observables have reached a 
plateau (as a function of the Monte Carlo time). 



4 Scaling and Density Analyses 

The analysis focuses on the odd sector of the RSIM, which is connected to the even sector 
through the standard scaling relations and their logarithmic counterparts [5j [6]. This 
connection is used to determine the specific-heat exponents through the density of zeros, 
in addition to a detailed analysis of the susceptibility and the Lee- Yang zeros. 

The FSS analyses of the magnetic susceptibility and Yang-Lee edge focuses on the 
correction-to-scaling exponents. From (12.41) . the reduced temperature is expresed in terms 
of the correlation length near criticality as 



£~ W (ln£ 



(4.1) 



Substituting (14. ip into (I2.3P and (12.61) . gives the scaling behaviour for susceptibility and 
the lowest lying Lee- Yang zeros in terms of the correlation length. Because there are no 



Table 3: The estimates for the leading-exponent ratios j/v and A/V from fits to the 
scaling behaviour of the susceptibility and first Lee- Yang zeros as well as estimates for 
(wy — yv)jv. These estimates agree with the theoretical values, which are 7/4, 15/8 and 
0, respectively. 



7/1/ 

(1/7 — yv)jv 

A/v 



p = 0.88889 



1.747 ±0.007 
-0.01 ±0.03 
1.879 ±0.004 



P 



0.75 



1.755 ±0.005 

0.02 ±0.03 

1.878 ±0.006 



P 



0.66661 



1.752 ±0.007 

0.01 ±0.03 

1.878 ±0.006 



logarithmic corrections to the FSS behaviour of the correlation length [6j, for sufficiently 
large lattices ^ may be replaced by L. This substitution then gives for the FSS of the 
susceptibility and the j th Lee- Yang zeros, 

Xl ~ L%(\nL) •> , (4.2) 

Tj ~ L~-(\nL)^^~ , (4.3) 

respectively. The values for the correction exponents 7 and v given in (12.81) and (12. lip 
have been numerically established for the RBIM in Refs. [T4"l [22] and that for and v in 
the RSIM has been verified in Ref. [27]. Therefore their confirmation in this setting as 
1/7 — yv = serves as a useful check on the accuracy of our method at the logarithmic 
level. 

The scaling analysis begins with p = 0.88889 (and /3 = 0.53781). A double logarithmic 
plot of xl against L is presented in Fig. la, and a fit to all data points yields 7/// = 
1.747(7) in agreement with (12.81) . The goodness of fit corresponds to a x 2 P er degree of 
freedom (% 2 /dof) of 0.8. The FSS correction exponent (1/7 — yi>) jv is extracted by a 
fit, corresponding to Fig. lb, of hix^ — 7/4 In L against In (InL), giving a slope —0.01(3) 
(with x 2 /dof ~ 0.8) compatible with the expected value of zero from (12.81) and (12.111) . 
These results are sumarized in Table [3] alongside the results of similar analyses at p = 0.75 
{j3 = 0.77125) and p = 0.66661 (j3 = 1.10). In each case, the x 2 /dof indicates a good fit 
and compatibility with the SSLJ theory is firmly established. 

Subleading corrections to (14.31) are expected to take the form 



Xl = AL^(lnL)" 



l + O 



1 



InL 



(4.4) 



The amplitude A of the leading term may firstly be estimated by fitting to xl — AL 1 ^. 
A subsequent fit to the parameters governing corrections to scaling in (I4.4p yields the 
estimate (1/7 - 7Z/) /u = -0.03(2) (with xV dof » °- 8 ) for at P = 0.88889. I.e., the inclu- 
sion of additive corrections does not lead to an improved estimate for the multiplicative 
logarithmic exponents. (We have also tested additive corrections of the form In InL/ InL 
in place of 1/lnL in ( 14.4D . While this leads to small improvement over ( 14.41) . it also does 
not significantly affect the estimates for the exponents of the multiplicative logarithms.) 




InL 



1.3 1.4 1.5 1.6 

ln(lnL) 



Figure 1: (a) FSS plot for xl at p c = 0.88889. The slope gives an estimate for 7/1/ of 
1.747(7). (b) Plot of Iv-Xl — 7/41nL against In (InL) giving slope —0.01(3), indicating no 
multiplicative logarithmic corrections to the FSS of the susceptibility. 



This observation holds for all quantities analysed below and for all values of p, and we 
henceforth refrain from reporting on additive corrections. 

The analysis of the leading FSS behaviour at p — 0.88889 of the first Lee- Yang zero 
of Ref. [3] is reconfirmed (to higher precision) in Fig. 2a, where, ignoring corrections, the 
slope of the log-log plot yields A/v = 1.879(4), compatible with the expected value of 15/8. 
The x 2 /dof here is 0.7. While the corresponding fit at p — 0.66661 yields A/v = 1.878(6), 
the equivalent fit (using all lattice sizes) for p = 0.75 gives the estimate 1.883(4), which is 
two standard deviations from the theoretical value. To avoid contamination due to scaling 
corrections, the smallest lattices may be removed from the analysis in standard fashion, 
and compatibility with theory is indeed restored when the smallest pair are removed from 
the fit; the L = 64 to L = 256 data again yields A/v = 1.878(6) (with x 2 /dof = 1.8). 
These reconfirmed results are also summarized in Table [3J 

To test the scaling relation (12.151) and the prediction (I2.16p . lnri + 15/8 InL is plotted 
against In (InL) in Fig. 2b for p = 0.88889. A fit to all points yields slope 



vA + At) 



v 



-0.01(2), 



(4.5) 



compatible with zero. The % 2 /dof here is 0.7. With the by now established theoretical 
values for A, v and v from (J2TB1 . fl2T9l) and fl2Tll . this yields 



A = -0.95(2) 



(4.6) 



a value compatible with (12.161) and confirming fl 2 . 1 5 [) in the RSIM. The analysis at p = 
0.66661 gives (z/A + Av)/v = —0.01(3), comparable to ( 14. 5ft and compatible with theory 
(the corresponding A value is —0.95(3)). At p = 0.75 one finds the estimate —0.04(2), 



: : 



u 




1.3 1.4 1.5 1.6 

ln(lnL) 



Figure 2: (a) FSS plot for the lowest Lee- Yang zero at p c = 0.88889. The slope gives an 
estimate for Ajv of 1.879(4). (b) Plot of lnri + 15/8 In L against ln(lnL) giving slope 
—0.01(2), indicating the absence of multiplicative logarithmic corrections to the FSS of 



using all lattice sizes. Again, compatibility with theory is restored by dropping the two 
smallest L values, where one again finds —0.01(3) (A = —0.95(3)). These estimates for 
A are summarized in Table HI 

Accepting, now, that the values A = 15/8 = 1.875 from ( 12.9ft and A = —15/16 = 
—0.9375 from (I2.16P are supported by the data for each dilution level, (12.101) and (I2.18P 
give the scaling exponents for the specific heat. In particular, the logarithmic correction 
exponent is identical to that of the density of zeros, 

a = a 2 ■ (4-7) 

Therefore a = 2 — 16a 2 /15 and a = a 2 can be extracted from the scaling form (12. 7p for 
the density of zeros, a form which is expected to hold for small enough r (i.e., close to 
the critical region). 

A robust method to determine the density of zeros from simulational data was devel- 
oped in Ref . [4j . Define the density of zeros along the singular line r > tyl {t) as 

g L (r)=L- d Y,Hr-r j (L)), (4.8) 

3 

where Tj{L) is the position of the j th zero for a lattice of extent L. Here j is called the 
index of the zero. Integrating along the locus of zeros gives the cumulative density of 
zeros (or index density) to be 



Gr(r) 



g L (s)ds = — for r^L) < r < r j+1 (L) 
o L a 



so that at a zero, it is given by the average 

GlML)) 



■2j 



2L d 



(4.9) 



(4.10) 




0.005 0.01 0.015 0.02 

rj(L) 



(b) 



o ■%, * * »<L# *^#f s 



N x * 



*x*$**¥^ ! 



a ® © *J® 



sifflxs stl, 



1.4 1.6 1.8 2 

In | lnr ■ (L) 



2.2 



Figure 3: (a) The integrated density of zeros at p = 0.88889 ('o' symbols), p = 0.75 
('x') and p = 0.66661 ('°'). The excellent data collapse (four data points for each of 
seven lattice sizes are involved at each value of p) and zero intercept indicate the correct 
transition point and the fitted curves give a 2 compatible with the expected value 16/15 
(a = 0). (b) The absence of multiplicative logarithmic corrections in the integrated 
density of zeros indicates that the specific heat logarithm exponent a = independent of 
the degree of dilution. 



From (12.71) . this may be fitted to the form 

G(r) = a 1 r a2 ()nrf 2 + a 3 . 



(4.11) 



allowing for an additional parameter 03 which determines the phase. A value of as greater 
or less than zero indicates that the system is in the broken or symmetric phase, respec- 
tively, so that a 3 = only at the transition point t = 0. The second criterion for a 
good fit is good data collapse. Note that this method does not allow an independent 
goodness-of-fit test [I]. 

The integrated density of zeros is plotted in Fig. 3a using the first four Lee- Yang zeros 
for lattices from size L = 32 to L = 256 (28 points in all) for each value of the dilution, 
demonstrating excellent data collapse in each case. For p = 0.88889, three-parameter fits 
to (14.111) . for small r, with a 2 = are supportive of the theoretical values of 02 = 16/15 
and 03 = 0; using the eight lowest data points, one obtains 03 = 0.000 000 2(4) and 
a 2 = 1.076(16). The latter results corresponds to the estimate a = —0.02(3), from (12.101) . 
The corresponding results in the p = 0.75 and p = 0.66661 cases are a 3 = —0.000 000 2(3), 
a 2 = 1.062(10) (a = 0.01(2)) and a 3 = -0.000 0001(4), a 2 = 1.066(15) {a = 0.00(3)), 
respectively. These results for a are gathered in Table HI 

We now accept that a 3 is indeed zero and the theoretical value a 2 = 16/15 (a = 0) 
from (I2.10p holds for each dilution. Potential multiplicative logarithmic corrections are 
detected by plotting lnG — 16/15 lnr against In (lnr) in Fig. 3b. A fit to all data points 
for p = 0.88889 gives a = a 2 = 0.012(3), which is four standard deviations from the 
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Table 4: The estimates for the specific heat exponent a and logarithmic-correction ex- 
ponents A and a. These are in agreement with the theoretical values from the strong 
universality hypothesis, namely a = 0, A = —15/16 and a = 0. 





p = 0.88889 


p = 0.75 


p = 0.66661 


A 
a 
a 


-0.95 ±0.02 
-0.02 ±0.03 
-0.02 ±0.05 


-0.95 ±0.03 

0.01 ±0.02 

-0.01 ±0.03 


-0.95 ±0.03 

0.00 ±0.03 

-0.04 ±0.05 



theoretical value of zero. However, focusing on the scaling region closer to the origin 
establishes compatibility with the theory. For example, fitting to the lowest eight data 
points yields a = d 2 = —0.02(5). The equivalent results for p = 0.75 and p = 0.66661 
are a = a 2 = —0.01(3) and —0.04(5), respectively. The corresponding fits are depicted 
in Fig. 3b and the estimates for a are summarized in Table HI These values constitute 
numerical evidence that a = a = 0, independent of dilution and in favour of strong 
universality. 

5 Conclusions 

The debate regarding the critical behaviour of the disordered Ising model in two di- 
mensions has persisted for over thirty years (most recently in (T8J [38]). Here the SSL J 
prediction [HI El [TO] for the multiplicative logarithmic corrections for the scaling behaviour 
of the susceptibility has been reconfirmed through a careful FSS analysis. In addition to 
this, the Lee- Yang zeros have been determined to high accuracy and their logarithmic 
corrections verified for the first time. 

The scaling behaviour of the specific heat in the site-diluted version of the model has 
been particularly difficult to pin down directly, and fits to the measured specific heat as 
a function of the double logarithm of the lattice extent [T5l [T8| [T9l [20] have been claimed 
not to be unambiguous |28j. Here an alternative approach has been taken, involving the 
density of Lee- Yang zeros. Using scaling relations [51 [6] to connect to the even sector 
of the model, the specific-heat scaling and correction exponents are clearly determined. 
Since the simulations are performed at three different values of the dilution (some quite 
large) the analyses presented herein for the susceptibility, the individual Lee- Yang zeros 
and for their densities, are unambiguously supportive of the strong scaling hypothesis. 
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